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Abstract.  A  simple  mesoscale  model  was  developed  for  discontinuous  dynamic  recrystallization. 
The  material  is  described  on  a  grain  scale  as  a  set  of  N  (variable)  spherical  grains.  Each  grain  is 
characterized  by  two  internal  variables:  its  diameter  and  dislocation  density  (assumed  homogeneous 
within  the  grain).  Each  grain  is  then  considered  in  turn  as  an  inclusion,  embedded  in  a  homogeneous 
equivalent  matrix,  the  properties  of  which  are  obtained  by  averaging  over  all  the  grains.  The  model 
includes:  (i)  a  grain  boundary  migration  equation  driving  the  evolution  of  grain  size  via  the  mobility 
of  grain  boundaries,  which  is  coupled  with  (ii)  a  dislocation-density  evolution  equation,  such  as  the 
Yoshie-Laasraoui-Jonas  or  Kocks-Mecking  relationship,  involving  strain  hardening  and  dynamic 
recovery,  and  (iii)  an  equation  governing  the  total  number  of  grains  in  the  system  due  to  the  nuclea- 
tion  of  new  grains.  The  model  can  be  used  to  predict  transient  and  steady-state  flow  stresses, 
recrystallized  fractions,  and  grain-size  distributions.  A  method  to  fit  the  model  coefficients  is  also 
described.  The  application  of  the  model  to  pure  Ni  is  presented. 

Introduction 

During  the  thermomechanical  processing  of  superalloys,  dynamic  recrystallization  controls  micro¬ 
structure  evolution  and  thus  important  aspects  of  the  mechanical  behavior  in  service.  For  instance, 
grain  boundary  migration  plays  an  important  role  because  it  is  one  of  the  main  phenomena  control¬ 
ling  the  final  grain  size.  The  present  paper  introduces  a  model  of  discontinuous  dynamic 
recrystallization  (DDRX)  which  describes  microstructure  evolution  for  a  distribution  of  isotropic, 
spherical  grains  of  various  diameters,  each  with  a  dislocation  density  p  which  is  the  same  for  all 
grains  of  identical  diameter  [1,2].  Topological  features  treated  elsewhere  [1]  are  not  considered  here 
in  order  to  simplify  the  approach. 

Description  of  Grain  Properties 

When  detenninistic  evolution  equations  ( i.e .,  with  no  stochastic  terms)  are  used,  all  grains  of  a  giv¬ 
en  age  r  have  undergone  identical  evolution  and  therefore  have  the  same  diameter  D  and 
dislocation  density  p.  Hence,  all  properties  of  the  grains  in  the  model  are  one-parameter  distribu¬ 
tions,  and  each  grain  is  characterized  by  its  age  r.  The  following  functions  can  then  be  introduced: 

-  the  number  of  grains  of  age  r,  N(r,t)  (number  per  unit  volume  and  age  time); 

-  the  plastic  strain  within  the  grain  s{r,t )  =  J  s(u) du,  in  which  the  strain  rate  £  is  assumed  to  be 

the  same  for  each  grain  (per  the  classical  Taylor  isostrain  crystal-plasticity  assumption); 

-  the  strain  hardening  of  the  grain  as  represented  by  its  dislocation  density  p{r,t )  (length  per  unit 
volume); 
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-  the  grain  diameter  D(  r,t). 

A  number  of  constraints  connect  the  various  functions;  e.g.,  the  overall  volume  is  constant  at  all 
times,  i.  e. , 

Vf,  |J'  N(t,  t)  D\r,t)  dr  =  1.  (1) 

Evolution  of  Grain-Property  Distributions 

Several  mechanisms  contribute  to  the  evolution  of  grain-property  distributions: 

(i)  Grain  boundary  migration.  A  simple  grain  boundary  migration  law,  based  on  differences  in 
strain  hardening  between  adjacent  grains,  or  more  specifically  between  a  specific  grain  and  the 
equivalent  homogeneous  matrix,  is  assumed,  i.e., 

(dD/8t)r  (r,t)  =  2 MT[_p{t)~ p{r,t)\,  (2) 

in  which  M  is  the  grain  boundary  mobility  and  T  is  the  line  energy  of  the  dislocations.  Such  a  fonn 
and  the  introduction  of  the  average  dislocation  density  p  are  justified  by  using  an  average  over  the 
matrix  surrounding  the  grain. 

Because  plastic  defonnation  occurs  without  volume  change,  p  must  be  a  surface  weighted  aver¬ 
age, 

p(t )  =  ^p(^T,t)N^T,t)D2(yT,t)  dr  j  Jo  A(r,t)D2(r,t)  dr.  (3) 

As  a  consequence  of  grain  size  evolution,  “old”  grains  (which  have  large  dislocation  densities) 
decrease  in  diameter  and  eventually  disappear.  They  are  replaced  by  new  ones  as  described  in  sec¬ 
tion  (iii)  below. 

(ii)  Strain  hardening.  Strain  hardening  and  dynamic  recovery  are  assumed  to  be  described  for 
each  grain  by  the  same  equation,  dp/ds  =  //(/?),  in  which  H ( p )  can  depend  on  temperature  as 

well  as  on  strain  rate  (in  the  same  way  for  each  grain). 

Several  forms  of  strain  hardening  and  dynamic  recovery  law  can  be  used,  e.g.,  the  Yoshie- 
Laasraoui-Jonas  law,  H(p)  =  h-rp;  the  Kocks-Mecking  law,  H( p)  =  h  - r^fp;  the  power  law, 

//(/?)  =  Hv+l  / pv .  In  each  case,  the  various  parameters  are  temperature  and  strain-rate  dependent. 

A  variant  exists  which  takes  into  account  the  additional  softening  due  to  grain  boundary  migra¬ 
tion;  but  this  law  is  then  no  longer  a  local  one  because  it  also  depends,  via  the  grain  boundary 
migration  velocity,  on  the  average  dislocation  density  [1-3], 

(iii)  Nucleation  of  new  grains.  The  nucleation  of  new  grains  (age  r  =  0)  is  postulated  to  be  pro¬ 
portional  to  the  grain  boundary  surface  S(t)  and  is  of  the  form  (dN/dt^  (0,/)  = 
f{yp)S{t),  where  t  is  the  current  time  and  S  is  the  grain  surface  area  per  unit  volume;  in  practice, 
this  latter  function  is  specified  as  /(/?)  =  kNpp.  The  exponent  of  the  Derby  relationship  between 
grain  size  and  stress  is  related  to  the  exponent  /;;  the  unique  integer  value  for  p  compatible  with 
experimental  Derby  exponent  is  3  [3].  In  the  present  work,  p  =  3  is  assumed,  although  p  could  be 
fitted  for  specific  sets  of  data. 

Using  the  above  equations,  the  change  with  time  of  the  grain-size  and  other  distributions  are  de¬ 
termined  as  shown  schematically  in  Fig.  1. 


2 


Figure  1.  Schematic  evolution  of  the  grain  size  distribution  in  the  DDRX  model 


The  flow  stress  a  associated  with  this  grain-size  distribution  is  assumed  to  be  given  by 
<j  =  ajub^Jp. 

Steady-State  Behavior 

When  the  material  is  strained  at  a  constant  strain  rate  and  temperature,  a  steady  state  is  obtained 
after  some  time  (or  strain)  during  which  new  grains  continuously  replace  old  ones. 

The  distribution  N{r,t)  is  then  constant  (i.e.,  independent  of  t),  and  the  diameter  distribution  is  a 

function  of  the  average  dislocation  density  according  to  a  specific  equation  (which  depends  only 
on  p,  the  dependence  on  s  and  temperature  being  implicitly  included  in  p).  Under  such  a  steady- 
state  condition,  Eq.  2  can  be  integrated  to  give: 

D(s,p)  =  2MT[_ps-F(s,p)\l£,  (4) 

in  which  F{s)  is  the  strain  hardening  and  dynamic  recovery  potential,  defined  by  p{s)  = 
(( dF(e,p)/de)_. 

The  steady-state  distribution  of  grain  size  g(D)  associated  with  such  a  one-parameter  distribu¬ 
tion  can  be  obtained  from  g(  D)  =  (dD/dt)  -{dDjdt')X_  j  rcnd,  in  which  rcnd  is  the  lifetime  of  the 

grains  and  the  indices  +  or  -  are  associated  with  growing  and  decreasing  grains,  respectively.  The 
above  distribution  necessarily  tends  to  infinity  when  D  — »  Dmax.  More  precisely,  a  g(D)  equivalent 

is  then  proportional  to  l/ ^Dm,a  -  D.  Nevertheless,  all  average  values  that  can  be  deduced  from 

Eq.  4  ( e.g .,  the  average  grain  size)  do  not  necessarily  have  the  same  problem  in  their  calculation. 
Moreover,  such  a  behavior  can  be  avoided  by  considering  that  the  maximum  diameter  is  likely  to  be 
not  precisely  defined  due  to  variability  in  the  material,  e.g.  variations  in  the  hardening  parame¬ 
ter  h  because  of  its  dependence  on  crystallographic  orientation. 

Determination  of  the  Model  Parameters 

For  given  strain  rate  and  temperature,  the  model  parameters  are  determined  from  stress-strain 
curves  measured  during  mechanical  tests,  such  as  hot  torsion  tests,  conducted  at  constant  strain  rate 
and  temperature.  After  a  pseudo-elastic  initial  domain,  the  initial  grains  strain  harden  before  the 

influence  of  the  newly  generated  grains  becomes  noticeable.  For  £e,  cre  =  apb^[p~()  where p0  is  the 


3 


initial  value  of  p.  The  first  part  of  the  stress-strain  curve  can  thus  be  used  to  fit  the  parameters  in¬ 
volved  in  the  definition  of  the  strain  hardening  and  recovery  function  //(/?);  see  (ii)  and  Fig.  2  for 
an  example. 


Figure  2.  True  stress-strain  curve  for  pure  nickel  determined  via  torsion  testing  at  900  °C  and 
0.1  s  1  and  model  fits  for  the  Yoshie-Laasraoui-Jonas  (YLJ)  and  power-law  (PW)  equations  using 
flow-stress  data  between  se  and  Sl 


In  the  second  step  of  the  fitting  procedure,  the  nucleation  parameter  kN  and  the  grain  boundary 
mobility  M  can  be  analytically  derived  from  the  steady  state  stress  crs  and  the  measured  average 
steady-state  grain  size  Ds  using  the  power  law  for  H ( p )  [3],  i.e., 

MT  =  ( 2v  +  3)(v  +  \)v+\v  +  2)  v  1  Hv+lsDs ( apb /  a s f{r+2\  (5) 


3(3v  +  4)(3v  +  5)(  v  + 1  )v+2Hv+l£ 


8(2v  +  3)2(v  +  2)l+1Ds2 


(apb/as) 


2(p+v+\) 


(6) 


In  the  case  of  other  formulations  for  H  ( p ),  these  values  can  still  be  used  for  M  and  ks  as  rough 

estimates  or  initial  guesses  for  performing  inverse-method  parameter  identification. 

Figures  3-5  present  model  predictions  of  the  evolution  (with  strain)  of  flow  stress,  average  grain 
size,  and  recrystallized  fraction  for  the  hot  defonnation  of  high-purity  nickel.  For  the  power-law 
(PW)  and  YLJ  curves,  the  values  of  M  and  kn  were  derived  from  Eqs  5  and  6  while  they  were  op¬ 
timized  for  the  “best-fit”  prediction.  The  fitted  curves  show  reasonable  agreement  with  the  available 
data  for  pure  nickel.  Furthennore,  the  value  of  the  exponent  p  in  the  nucleation  equation  can  be 
tuned  in  the  neighborhood  of  p  =  3  to  improve  the  agreement. 

The  model  can  also  be  used  for  more  detailed  predictions  such  as  the  grain  size  distribution,  as 
illustrated  in  Fig.  6.  As  mentioned  previously,  the  exact  shape  of  this  distribution  can  vary  due  to 
variations  in  the  material  parameters;  the  inclusion  of  local  neighborhood  effects  (which  would  re¬ 
quire  a  more  complex  model  formulation)  may  also  lead  to  different  fonns  of  the  grain-size 
distribution  that  show  more  of  resemblance  to  experimental  observations. 
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Figure  3.  Comparison  of  predicted  and  measured  flow  stress  of  Ni  deformed  at  900  °C,  0.1  s  1 


Figure  4.  Predicted  vs  measured  average  (volume  weighted)  grain  size  of  Ni  (900°C,  0.1  s  '  ) 
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Figure  5.  Predicted  recrystallized  fraction  for  pure  nickel  deformed  at  900  °C  and  0. 1  s~ 1 
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Figure  6.  Predicted  (YLJ  equation)  steady-state  grain-size  distribution  for  pure  nickel  defonned 
at  900  °C  and  0.1  s”1 


Closing  Remarks 

The  DDRX  model  can  also  be  used  for  predicting  the  behavior  of  alloys  during  hot  working.  When 
sufficient  data  are  available,  it  can  be  used  in  the  same  way  as  above  for  pure  nickel.  In  the  case  of 
solid  solution  nickel-niobium  alloys,  for  example,  analytical  descriptions  of  the  strain  hardening 
and  dynamic  recovery  parameters  have  been  established  [4].  The  grain  boundary  mobility  and  pos¬ 
sibly  the  nucleation  parameter  are  affected  by  solutes  and  second  phase  particles.  The  former  effect 
has  been  addressed  in  a  companion  paper  in  this  proceedings  volume  [5]. 
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